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• accretion disk during the radio plateaux is always associated with radiation- 

driven wind. The internal shock forms in the previously generated slowly 
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1 Introduction 



In active galactic nuclei (AGN) & quasars, large superluminal outflows/jets are known to 
exist for almost half century yet remain poorly understood. This is probably due to the 
lower characteristic accretion disk temperature as well as the fact that the central part of 
these systems is usually obscured by a large amount of dust. The long time scales associated 
with these massive systems pause an additional observational problem. Discovery of the first 
microquasar GRS 1915+105 in 1994 in our Galaxy revived great hope to study the disk-jet 
connection in general and the origin of superluminal jets in particular [Ij. Microquasars are 
closer, smaller and show faster variability that is easily observable. GRS 1915+105 has shown 
exceptionally high variability in X-rays as well as in radio. Since 1996, rich X-ray variability of 
this source is observed by RXTE [2ii3j and by the Indian X-ray Astronomy Experiment, IX AE 
[4]. Belloni et al. (2000) have classified the complex X-ray variability of GRS 1915+105 in 12 
separate classes on the basis of their light curves and color-color diagrams and suggested three 
basic states of this source, namely hard state and two softer states with different temperatures 
of the accretion disk. 

The radio emission from GRS 1915+105 can be broadly put into two classes; (1) radio 
emission close to the compact object (j200 AU), and (2) radio emission at large distances (> 
240 AU). The former class includes (a) steady radio jets (radio plateau state), (b) preplateau 
flares, and (c) oscillations/baby jets (discrete jets) of 20-40 min duration in infrared (IR) & 
radio while the latter class includes large superluminal radio jets. The steady radio jets of 20 
- 160 mJy flux density are associated with the canonical low hard X-ray state and observed 
for extended durations [6|,|7j. These are optically thick compact jets with velocity f3 of 0.1-0.4 
P, [9j. The radio emission is correlated with the X-ray emission as Lradio oc L^^ for several 
different sources [10]. Pooley & Fender (1997) observed radio oscillations with delayed emission 
at lower frequency. Simultaneous X-ray, IR and radio multi-wavelength observations provided 
first major step in our understanding of disk-jet interaction and suggested that the spike in 
X-ray coincides with the beginning of IR flare and it has self absorbed synchrotron emission 
associated with adiabatic expansion [T^lIISl- These are also compact jets with velocity /3 ~ 1 
[S]. The relativistic superluminal jets with up to 1 Jy flux density have steep radio spectrum 
and are observed at large distances few hundred AU to 5000 AU from the core lU [HI |8] . These 
radio jets are very energetic with luminosity close to the Eddington luminosity, LEdd and have 
been observed in several sources [HI [151 HSl [HI- Progress in our understanding of these jets, 
especially of their connection to the accretion disk has been slow. The physical connection 
between X-ray emission and the superluminal flares has been the hardest to understand. 

In this paper, we investigate the association of large superluminal jets with the radio 
plateaux. We have analysed the available RXTE PCA/HEXTE X-ray data during radio 
plateaux and the radio flare data from the Green Bank Interferometer (GBI). Yadav (2006) 
has described the analysis procedure and the selection criteria in detail. 
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a: Integrated 3-150 keV X-ray flux in 10 ergs cm & 



b: VLBA radio data (see text for details) 

Table 1: All selected superluminal radio flares and their properties (see text for details). X-ray 
Properties from RXTE PCA/HEXTE data during preceding radio plateau. 

2 Observations and analysis 

The typical sequence of events for a superluminal radio flare is shown in Figure [1] for the 550 
mJy radio flare on 1997 October 30 (gap in GBI data at peak) and for the 340 mJy radio 
flare on 1999 January 30. The start of the radio flares is offset to zero. A superluminal event 
starts with small preplateau flares followed by a steady long plateau, followed by superluminal 
radio flares. The preplateau flares are discrete ejections of adiabatically expanding synchrotron 
clouds with flat radio emission and are similar to the oscillations/discrete jets. The exponential 
decay is an important characteristic of superluminal radio flares which differentiates them from 
the other class of radio flares which occur close to the compact object [20]. It is also found 
that the radio plateau is always associated with a superluminal radio flare [211 122] . We 
searched 2.25 GHz GBI radio monitoring data during the period from 1996 December to 2000 
April and selected radio plateaux and the following radio flares which decay exponentially. All 
the selected superluminal flares are given in Table [T] along with associated X-ray properties. 
One more radio flare on 2001 July 16 which was observed by the Very Large Baseline Array 
(VLBA) and Ratan radio telescope is also added in Table [H VLBA observations clearly 
showed an ejecta well separated from the core [23] . 

We study X-ray properties during radio plateaux within the preceding week from the start 
of superluminal flares to avoid changes in accretion disk over the long durations of plateaux 
(from -7 to -2 in Figured]). We also avoid the last day of plateaux as radio data suggest rapid 
changes in the accretion disk. For the timing analysis, we used single-bit-mode RXTE/PCA 
data in the energy ranges 3.6 - 5.7 and 5.7 - 14.8 keV. Normalised power density spectra 
of 256 bin are generated and co-added for every 16 s. During the plateaux, it is in a very 
high luminosity state (VHS) with power law index P > 2 . The spectrum is dominated by 
the Compton scattered emission (> 85%) rather than the disk. A model consisting of disk -|- 
power law -|- a comptonised component (CompTT) with a gaussian line at 6.4 keV is used for 
the X-ray analysis. The integrated radio flux is calculated by integrating the fltted exponential 
function over a duration three times the decay time constant. Further details of X-ray and 
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Figure 1: GBI 2.25 GHz radio data for the 1999 January 30 (start at MJD 51204.7) and 1997 
October 30 (start at MJD 50750.6) large superluminal radio flares. Start of the flares is offset 
to 0. The dotted line connects the data points of the 1999 January 30 superluminal radio flare 
to guide the eyes. 



radio data analysis are given in Yadav (2006) [20] . 



3 Results and discussion 

It is increasingly believed that the coronal material (and not the disk material) is ejected 
prior to radio flares [2H |22l [181 ES]. It is consistent with the observation that radio flares 
are observed during the transition from the low hard state to the high soft state but never 
observed during the opposite transition [21]. It also agrees with the suggestion that the spike 
in X-ray during the (3- class is associated with the change of the X-ray state due to a major 
ejection episode [2S]- Once this coronal material is ejected, it invariably decouples from the 
disk. However, this decoupling may not produce any observable effect during the plateaux 
as inflow and outflow in the accretion disk are in equilibrium. Both the accretion disk and 
the radio flare are in steady state. On the other hand, decoupling is supposed to produce 
maximum effect during the superluminal radio flares which are observed at large distances of 
a few hundred to few thousand AU from the core [6l[Tll[8]. This has made disk-jet connection 
the hardest to understand for superluminal radio jets [IHl [I9j. 

In Figure [21 we plot the total 3-150 keV X-ray flux during the preceding plateau vs the 
integrated radio flux of the superluminal flare. We derive a correlation coefficient of 0.99 sug- 
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Figure 2: Correlation between the total 3-150 keV X-ray flux from spectral analysis of RXTE 
PCA/HEXTE X-ray data during preflare plateau state and the integrated flux of the super- 
luminal radio flare. The solid line is the linear flt to the data points. 



gesting a strong connection between the total X-ray flux during the preceding plateau state 
and the integrated radio flux of the superluminal flare. The total X-ray flux is dominated 
by the Compton scattered emission or coronal emission (> 85%). Using the total Compton 
scattered emission flux instead of the total X-ray flux improves this correlation with a correla- 
tion coefficient of 0.996. In Figure [3], we plot the QPO frequency from our timing analysis of 
X-ray data during the plateau state as a function of the decay time constant of the following 
superluminal radio flare. This also shows a tight correlation, with correlation coefficient of 
0.98. The QPOs are believed to be associated with the coronal flow. The remarkable feature 
of our flndings here is that the parameters calculated using completely independent spectral 
and timing analysis bring out a clear connection between the accretion disk during the plateau 
state and the following superluminal radio flare. 

In the last column of Tabled! we give the calculated absorption column density Nh- The 
Nh ranges from 10x10^^ to 15x10^^ cm~^ which are higher than the commonly used Nh ~ 
5 - 6x10^^ cm~^ for spectral analysis of GRS 1915+105 . The calculated Nh shows a tight 
correlation with the CompTT flux and not with the disk blackbody flux. This rules out the 
possibility that the enhanced Njy may be due to overestimation of the disk normalisation. 
Yadav (2006) has discussed this in detail and compared these with results obtained from 
Chandra and ASCA data [271 [28]. Lee et al. (2002) have analysed Chandra & RXTE X- 
ray data during a radio plateau state and have suggested presence of disk wind. Kotani et 
al. (2000) also came to a similar conclusion using ASCA data. The wind terminal velocity is 
estimated to be of the order of 10*^ cm s~^ {(3 < 0.01). The lower limit of bolometric luminosity 
Lhoi ~ Lx = 6.4x10^^ ergs s~^ which is 0.35 of Eddington Luminosity LiEdd, for a black hole 
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Figure 3: Correlation between the quasi-periodic oscillation (QPO) frequency from timing 
analysis of RXTE/PCA X-ray data during prefiare plateau state and the decay time constant 
of the superluminal radio flare. The dotted line is the linear fit to the data. 



of mass 14 p7] . The corresponding lower limit of bolometric luminosity for the X-ray flux 
listed in Table [1] falls in the range 0.35 - 0.48 of hEdd- When GRS 1915+105 is accreting near 
IjEdd^ the presence of a radiation-driven wind is always expected and the wind density should 
be a strong function of the disk luminosity. Our derived values of show strong dependence 
on observed total X-ray flux with correlation coefficient of 0.995 |30j . 

The calculated power density spectra shown in Figure H] also independently lend support 
to the presence of wind during the plateau state. Figure H] shows PDS spectra observed on 
2001 July 11 and 1998 April 28. The Compton scattered flux increases from 1.5x10"^ ergs 
cm~^ s~^ on 2001 July 11 to 1.9x10"^ ergs cm~^ s~^ on 1998 April 28 while the observed 
QPO frequency decreases from 2.4 Hz on 2001 July 11 to 1.4 Hz on 1998 April 28. It is clear 
from Figure m that the power at frequency v > 0.2 Hz is less in the PDS observed on 1998 
April 28 than that observed on 2001 July 11. The fast variability is suppressed by photon 
scattering in the enhanced wind on 1998 April 28, hence reducing the power in the PDS at 
higher frequencies. Shaposhnikov & Titarchuk (2006) have discussed the decrease in the PDS 
power observed in Cyg X-1 at higher frequency [u > 0.1 Hz) as the wind increases. In Cyg X- 
3, the PDS power in the low hard state drops to below 10^^ (rms/mean)^ Hz~^ at frequencies 
u > 0.1 Hz as a dense wind from the companion always envelops the compact object [50] . 

These results support the internal shock model for the origin of superluminal flares |31j . 
The internal shock should form in the previously generated slowly moving wind from the 
accretion disk with (3 < 0.01 as the fast moving discrete jet with /? ~ 1 [8j catches up and 
interacts with it. Both the components, slow moving wind and fast moving jet, are related 
to the accretion disk during plateau state and the strength & speed of these two components 
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Figure 4: Normalised power density spectra in 0.09 - 9 Hz frequency range observed on 1998 
April 28 and 2001 July 11. Strong QPOs are seen at 1.4 Hz and 2.4 Hz in PDS observed on 
1998 April 28 and 2001 July 11 respectively. The total Compton scattered flux is also given. 



will determine the power of the internal shock. The wind deposits a large amount of energy 
as rhwind approaches Thaccr prior to the switch-on of superluminal flare [20]. Thus, the internal 
shock model can easily accommodate high jet power requirement > 10^^ ergs s~^ [13] and can 
explain the shifting from thick to thin radio emission during superluminal flares [31]. Our 
results in Figure [5] which shows the peak flux of superluminal radio flares as a function of 
strongly support our description of superluminal radio flares. A fit to the data (dotted line) 
suggests that for wind strength corresponding to N/^ < 8.3±1.5xl0^^ cm~^, no superluminal 
jet will be produced. The absence of superluminal jets during the class /3 in GRS 1915+105 is 
attributed to the absence of wind [20]. The calculated during class P is below this critical 
Nh value [30]. This model of superluminal jets can provide simple explanations for (1) why 
the direction of superluminal jets may differ from the direction of the compact jets, (2) why 
the phase lag is complex during plateaux, and (3) why compact (discrete) jets are absent in 
Cyg X-3 [20J. 

In Figure [6], we show a sketch of the heliosphere due to disk wind and mark the locations 
of various types of radio emission seen in GRS 1915+105 . All the radio flares observed in 
GRS 1915+105 can be broadly put into two groups on the basis of their flux, radio spectrum 
and spatial distribution; (1) the superluminal flares (200-1000 mJy) which have steep radio 
spectra and are seen at large distances (> 240 AU), and (2) all other flares (5-360 mJy) 
which include the preplateau flares, radio oscillations & discrete flares and the steady radio 
emission during the plateaux. All these flares have flat radio spectra and are observed close to 
the compact object. All these radio flares are consistent with the ejection of an adiabatically 
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Figure 5: Correlation between the absorption column density, Nh and the peak flux of the 
superluminal radio flares. The peak flux of a superluminal flare is calculated using exponential 
profile fitting if there is a gap in radio data. The dotted fine is a linear fit to the data. 



expanding self absorbing synchrotron cloud from the accretion disk and this cloud is supposed 
to consist of coronal mass. During the radio plateau, it may be confined expansion due to the 
presence of a dense wind which agrees with the source size ^8j. The IR & radio oscillations 
(5-150) are periodic ejections of adiabatically expanding self absorbing synchrotron clouds 
[T3t [32] . The preplateau radio fiares (50-360 mJy) are discrete ejections which are closely 
spaced in time and hence produce overlapped radio fiares. The preplateau fiares have been 
modeled as adiabatically expanding self absorbing clouds ejected from the accretion disk (like 
in the case of oscillations & discrete jets) which explain reasonably well all the available 
observed data of time delay for radio emission at lower frequency All these types of radio 
emission in GRS 1915+105 starts somewhere close to the corona. 

As discussed above, the rhaccr, as inferred from the X-ray fiux, is very high during the radio 
plateaux and Lboi approaches LEdd- It is suggested that such hot accretion disk during the 
radio plateaux always accompanies with a radiation-driven wind. The Nj^^ is tightly correlated 
with the COMPTT fiux. This is analogous to the solar wind originating from the dense solar 
corona. Since the disk wind power & speed are higher than those of solar wind (the average 
solar wind speed ~ 450 km/s), the size of heliosphere due to the disk wind may be around 
200-300 AU (the solar heliosphere size is of 100-150 AU and it is supposed to vary with the 
solar activity). It is expected that the size of disk heliosphere will increase as the wind power 
increases with disk luminosity as discussed earlier. As shown in Figure El superluminal radio 
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Figure 6: Sketch diagram (not to the scale) of the hehosphere created by disk wind similar 
to the solar heliosphere. The locations are marked where different types of radio jets are 
produced. Superluminal jets start at the boundary of the disk heliosphere while all other 
types of the jets start somewhere close to the corona. 



jets are suggested to appear at the boundary of disk heliosphere. It is expected that strong 
superluminal jets should appear at larger distances than the distance where weak superluminal 
jets appear. A superluminal jet with peak flux of ~ 200 mJy was reported around 240 AU on 
2001 July 16[23] while strong superluminal jets are suggested to appear around 500 AU (the 
peaks of these superluminal flares are not observed but the data are consistent with the above 
suggestion) [8] . This can be easily tested in future with VLBA and other high resolution radio 
data if we catch the start of superluminal flares in GRS 1915+105 as well as in other LMXBs. 

4 Conclusions 

We have provided a tight correlation between accretion disk and superluminal jet parameters. 
We find that rriaccT during the plateaux is very high and suggest that the accretion disk during 
the radio plateaux is always associated with a radiation-driven wind. The internal shock forms 
in the previously generated slowly moving wind (during plateau) with /3 < 0.01 as the fast 
moving discrete jet (usually at the end of plateau) with /5 ~ 1 catches up and interacts with 



9 



it. The strength & speed of these two components determine the power of the internal shock 
as well as of the superluminal jets. The peak flux of the superluminal jets is tightly correlated 
with the wind power. We discuss the implications of this work and some of these can be 
checked in future. 
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